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effect of this difference is included in the relativistic correction. We attempt 
to define the difference in terms of NRQCD matrix elements and to sepa- 
rate the effect of the difference in several decays. It turns out that the total 
relativistic correcttion can be estimated by the difference. Numerically the 
correction will enhance the decay widths considered here. 
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Quarkoniuin systems generally are thought as simpler than light hadrons and it may be 
easier to handle them in the framework of QCD. However, only recently we have been able to 
treat their decays and productions rigorously, based on a factorization with non-relativistic 
QCD(NRQCD) [|^], where the effect of short distance is handled with perturbative QCD and 
the effect of long distance is parameterized with NRQCD matrix elements. The factorization 
is performed by utilizing the fact that a heavy quark moves with a small velocity v in the 
quarkonium rest frame and an expansion in v can be employed. In this factorization, a 
quarkonium decay, e.g., like J/\l/, can be imagined at the leading order of v as the following: 
The c- and c- quark in J/\l/ has certain probability to be freed at the same space point and 
this cc pair decays subsequently. The probability has a nonperturbative nature and is at 
order of f °, while the decay of the cc pair can be treated with perturbative QCD. Various 
effects at higher order off can be taken into account, e.g., relativistic effect, and the effect of 
that the freed cc pair does not possess the same quantum numbers as those of J /'^ . In order 
to make the factorization consistently the quarkonium mass is forced to be approximated 
as twice of the heavy quark mass. In general the quarkonium mass is not the same as twice 
of the heavy quark mass, the difference, we call it as the mass difference, is at order of v"^ 
and is made by nonperturbative physics, which binds the heavy-, antiheavy quark and other 
possible light dynamical freedoms into a bound state. This difference makes the phase space 
of a quarkonium decay differently than that of a heavy quark pair decay and introduces a 
correction at order of v"^. The situation here is relatively similar as this in deeply inelastic 
scattering, where the mass of the initial hadron introduces the target mass effect . This 
effect can be thought as a kinematical effect at higher- twist, the other higher twist effects 
are dynamical, which come from the fact that partons inside a target can have a transverse 
momentum and can be off-shell, ■ ■ ■ etc., these together are corrections to the parton model. 
In a quarkonium decay, like J/\l/-decay, the correction from the next-to-leading order of v 
is characterized by a single matrix element in NRQCD and both kinematic- and dynamical 
effects are included in the correction, as shown below explicitly in considered processes. 
However, the kinematic effect can be separated from the high-order corrections, and it can 
be interpreted in term of the mass difference. A decay width obtained with and without the 
separation should be at the same accuracy at the considered orders. This fact leads to that 
the effect at the next-leading order of v"^ can be interpreted totally by the mass difference 
in the decays considered here. In this work we will first try to define the mass difference in 
terms of NRQCD matrix elements and then calculate the correction from the next-to-leading 
order with and without the separation. 

We start with the energy momentum tensor T^y of QCD, a nice discussion about proper- 
ties of the tensor and interesting references can be found in 0. For our purpose we consider 
the full QCD containing gluons, Nf flavors of light quarks and one flavor of heavy quark Q. 
The tensor then can be written as 

T^'^{x)=T^\x)+T^\x)+T!;^{x) (1) 

where Tq" is the contribution from the heavy quark, Tq^ is from gluons and Tj^'^ is from 
light flavors. For a state of a hadron with momentum p one has by definition: 

{p\jd'xT^^{x)\p) _ 



If the hadron is at rest, the zeroth component in the right side of Eq.(2) is the mass M 
of the hadron. The energy momentum tensor can be decomposed into a trace part and a 
traceless part, 

T^^(a;) =f^'^(a;)+t^^(x), (3) 

(7^,T^-(x) = 0. (4) 

From the space-time covariance one can derive 

P^{p\T^^u\p) =P,,Pu, (5) 

p'{p\f,,\p)=p^p,-^g^M'', (6) 

p'{p\f,M = \9,uM' (7) 

where the factor p^ appears because we take the normahzation of a state as {p\p') = 
{2tt)^6^{p — p'). To interpret the mass difference we match the heavy quark fields of full 
QCD into NRQCD fields. We perform this matching at the tree-level and for consistency 
we neglect the contributions from the trace-anomaly and the anomalous dimension of quark 
masses. The heavy quark contribution Tq^ in the full QCD reads: 

T^ = Tl^" + ^\ (8) 

f^" = \g'''mQQQ. (10) 

In the above Q is the dirac field for the heavy quark Q, D'^ is the covariant derivative. The 
last term with g^"^ in Tq^ can be simplified with the equation of motion. We take the state 
\p) as a quarkonium state \H) and match the matrix element of Tq" into NRQCD matrix 
elements, we obtain: 

Mh = 2mQ - {H\Tk\H) + ^ mf{H\qjqj\H) + 0{v^) + 0(a,) (11) 

/ 

where T^ is an operator defined in NRQCD, ipq and xq ^"^^ the quark and antiquark fields 
in NRQCD respectively, Mq is the pole mass of the heavy quark. The operator Tk measures 
the sum of kinetic energies of a heavy and antiheavy quark. For H with a nonzero spin the 
average of the spin is implied. The summation is over all light fiavors. The matrix elements 
{H\qfqf\H) are at least at the order of f ^. This can be estimated by taking a QQ pair instead 
of the state H, because the operator gjgj is a color-singlet, the matrix element is nonzero 
provided that at least two gluons must be exchanged between the operator and the QQ state. 
A gluon coupled to Q or Q gives a factor v. The contribution from this operator to M^ 
is suppressed by a factor rrif/mQ, relatively to the contribution from Tk, can be neglected 



^r = niQi^D^^" + ^D^i^Q - T^g'^'Q^i ■ DQ}, (9) 



for charinonium as the mass ratio is small enough. For bottonium system, because charm 
quarks are included in the full QCD, neglecting this contribution from charm quark may be 
problematic, because the mass ratio is not so small. However the matrix element {H\cc\H), 
where c and c is the dirac field for charm quark, is suppressed by w^AQf.Q/m^ at least, hence 
the contribution from this operator suppressed by AQQ-^/mbrric relatively to that from Tk- 
Giving the fact that Aqcd as a characteristic scale for long-distance physics is at order of 
several hundred MeV, this contribution can be neglected too. 

Neglecting the contributions proportional to rrif and other higher order effects we obtain 

MH = 2mQ-{H\TK\H). (13) 

Taking this result we obtain another relation in the quarkonium rest frame within our 
approximation: 

2{H\Tk\H) = -{H\T^ + T°°|i7). (14) 

This relation may be regarded as a field theory version of the virial theorem of quantum 
mechanics, which relates the expectation value of the kinetic energy to the expectation value 
of the potential energy. The above equations hold only at the tree-level approximation in 
field theory. If one takes the contributions from the trace-anomaly, the contributions from 
the anomalous dimension of quark masses and those from the matching beyond the tree- 
level, into account, these relations become more complicated. For higher exited states of 
quarkonia the relation needs to be corrected by higher orders of f ^ as v"^ becomes larger for 
these states. The matrix element {H\Tk\H) is positive and it indicates that Mh < 2mQ. 

If the mass of a quakonium and the pole mass of the heavy quark are known, one may 
determines how large the matrix element {H\Tx\H) or the binding energy is. We take T and 
J/\l/ as examples. The pole mass of b-quark is determined precisely from a study of lattice 
QCD [^, whose result is rrih = S.OGeV. For charm quark we take the pole mass determined 
from the D meson semileptonic decay @], the value is rric = 1.65GeV. We obtain: 

{T\Tk |T) ^ 0.54GeV, (15) 

(J/^|Tx| J/^) ^ 0.20GeV. (16) 

With these values one can also determine the velocity of b- or c-quark in T or J/\l/ respec- 
tively, by identifying {H\Tk\H) = rnqv'^. We obtain f^ ^ 0.1 for b-quark and f^ ^ 0.12 for 
c-quark. The velocity for c-quark is smaller than that from usual estimation. For c-quark it 
should be noted that a precise determination of the pole mass is still lacking and the higher 
order effect in Eq.(13) may be significant. 

Now we consider the decay J/\E' -^ d^d" . Starting with the relevant ^-operator we 
obtain the ^'-matrix element 

{t{ki)A-{k2)\S\J/^{P)) = -te^Q,{2nf6\P-h-k2) 

-0y,u{k,hAk2)^2^r{YT{q,P)) 

Ti^{q,P) = Jd'xe"^--{0\c{x),c{x)i\J/^) (17) 



To perform the NRQCD factorization we match the nonperturbative object r{q,P) into 
NRQCD matrix elements by the small velocity expansion. Up to the order of f ^ it reads: 

r,,{q,P) = {27cy6\q-po)h^{P^l'P-UO\xW^c\Jm 

- A(2vr)^5^(g -po)(P+V^-)., — (0|(D^Xc)VV,| J/vI/) (18) 

where pg = {2mc, 0, 0, 0), P± = (1 ± 7°)/2. The first term is at order of v^, the remaining 
terms are at order of v"^. The relativistic correction is not only from these terms but also 
from P^ = (Mj/ij,, 0, 0, 0) in the S'-matrix element. For the terms at v"^ in r{q,P) it can 
be replaced by P = po because we neglect the effect at higher orders. For the first term in 
r(g, P) we can write the S'-matrix element with translational invariance as: 

{t{h),r{k2)\S\J/^{P)) = ie'Q,-^u{hWv{k2){27i)'6\P -k,- A;2)(0|xlaVe| J/^) 

+ 

+ ' (19) 

where stands for the contributions from the terms at f^ in T{q, P). The x-dependence 

in the NRQCD matrix element can be expanded. This expansion will lead to a tower of 
operators with n-derivative with x^. Similar cases in quarkonium productions are considered 
0, where one can add the effect due to these operators to cross-sections at the leading order 
of V and it results in that the hadronic phase-space is recovered. Here we start with the 
hadronic phase-space and try to approximate it with the partonic-phase space. After the 
expansion we obtain the S'-matrix element up to order of f ^: 

(£+(fci),r(A;2)|S|J/M/(P)) = te^QMhhHk^) 

■ {{2nr6\po - fci - k2)i-^{0\xy^Mm 

Atti^ 

+ i^^(2^)''^'(Po ~h- A;2)(0|(D2xc)VV,| J/vl/)} (20) 
With the above S'-matrix element we obtain the decay width 
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+ ^{{Jm^PlaXc\0) ■ (0|(D\,)V^,| J/vl/) + h.c.)] + 0{v') (21) 

where the spin average for J/\l/ is imphcit in the products of the matrix elements. This 
result agrees with that derived by the matching procedure in [|1|. The relativistic correction 
is represented by the product of matrix elements: 

a = {j/^\ijlaxc\0) ■ m-D\,yaij,\J/^) + h.c. (22) 

The decay width depends on Mj/^ only implicitly through the matrix elements, while origi- 
nally there is an explicit dependence on Mj/-^ in the term detailed in Eq.(19). However this 
term can be treated by the direct expansion of the ^-function 

SiMj/^ -kl- A;°) = 5{2m, - A;° - A;°) - ^Mj/^^5{2m, - fc? - A;°) + 0{v^) 

oki 

AMj/^ = Mjj^ - 2mc (23) 

where AMj/^ is at order of v"^ as shown before. With this expansion and r(g,P) given in 
Eq.(18) the decay width is 



,27r 1_ . , o,, AM.7 

mi rric 



T{JI^ ^ e*r) = a=Q?4^-i,||{0|xy*.| J/*>p(i + ii=^) 



+ ^«l + 0(.*). (24) 

Both results should be equivalent up to order of f ^, one can hence find a relation between 
a and Aj/^ and rewrites the decay width as 



.2vr 1 ,,, , X , , ,,,, 4AMj 



r( J/* ^ l!*r) = a^Q^,i^^\mxlaMJ/n'{l - :=E^) + 0{v'). (25) 



As discussed before, the quantity Mj/i^, is negative, hence the relativistic correction seems 
to enhance the decay width. This is in contrast to the model calculation where one uses 
a wave-function computed in nonrelativistic potential models. With the wave-function the 
parameter a can be estimated as 

a = ^R*{r)s7^R{r)\r=o + h.c. (26) 

where R{r) is the radial wave function. From above a is negative, the decay width therefor 
is reduced. But as discussed in ||l|, this parameter estimated in this way is not well defined 
because it is divergent if the potential is Coulombic at short distances. Certain subtraction 
is need to make it finite. It is unclear how to do the subtraction in a model calculation 
to make the estimation meaningful for the renormalized matrix element in a. In model 
calculations certain regularization is used to obtain a finite relativistic correction, like in [|^ . 
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In rewriting the decay width we used the relation between a and AMj/ij,. It should be 
kept in mind that this relation is not proved in general and it applies only in the detailed 
cases considered here. Similar analysis can be done for t],. — > 77 and the result is that the 
kinematic effect does not appear at the next-to-leading order of v. 

The other process to be considered in this work is J/\l' -^ light hadrons. It should be 
noted that for processes involving strong interaction NRQCD factorization is in general per- 
formed for squares of amplitudes. However, the vacuum saturation is a good approximation 
in NRQCD, this enables us to work at the level of amplitude. With the vacuum saturation 
and at leading order of Ug the process to be considered is J/^ -^ 3 Gluons. The S'-matrix 
element can be written: 

(3G|5| J/M/(P)) = z{2ttY5\P -h-h- h) J ^Jr{A{h, h, h, q)r\q, P)}, (27) 

rf^.(g) = Jd'xe''^-^0\c,{x)cm\Jm, (28) 

A{h,k2,h,q) = -^^fT«^T"^T«3^ i ^ 

(fci - q)^ - m^) (fci + fc2 - q)^ - m^) 

■ 7 ■ ^i'(7 ■ (^1 -q)+ m^h ■ e^Kl ■ ih + k2-q)+ m^)-f ■ ^^ 

+ Permutation of 1, 2, 3. (29) 

The indices 1,2,3 are labels for the 3 gluons. To calculate the decay width up to v"^ we 
match the nonperturbative object T^{q) into NRQCD matrix element and only keep the 
contributions up to order of v"^ and take the contribution at v"^ from the (5-function for the 
energy conservation into account. In the calculation of the effect at v'^ one should use physical 
spin-sums for gluons and be careful with the integration over the 3-body phase-space. The 
detailed calculation is tedious but straightforward. We obtain: 

20 
r(J/M/ ^ light hadrons) = ——aUmxlcr^c\Jm\\n^ - 9) 

+ (l^,r2-I).i-a} + 0(t;4). (30) 

492 2^ m2 ^ ^ ^ ^ ^ 

Similarly as for the leptonic decay we can directly expand the 5-function for the energy 
conservation and obtain the decay width as 

r{J/^ -. light hadrons) = ^,aUmxla^.\Jm\^' - 9 + ^^^3/^"'^^'^^ ) 
Both widths are correct up to order of v'^, we can hence rewrite the width as: 



r{J/^ -. light hadrons) = ^^,amxlcrA\J/n'[^' - 9 - (^tt^ - 7)^^] 

+ 0{v^). (32) 

If we take rric = 1.65GeV as before, the width for the decay into leptons and into light 
hadrons is enhanced at the level of 16% and of 48% respectively.. The hadronic width receives 
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a substantial correction. From Eq.(32) and Eq.(25) one can also obtain the corresponding 
widths for T. The same level of the enhancement is also found for T —>■ i'^i~ and for 
T — >■ light hadrons. To compare with experimental data we build the ratio: 

r(T ^ light hadrons) 
^^^ r(T -^ i+i') 

= 5040.0a^(mfe)(1.276 - 0.081a,(mb)) (33) 

where we have taken Mj- = 9.460GeV and the pole mass rrif, = S.OGeV. We have also added 
the one-loop QCD correction in each decay width p|Q. In this ratio the one- loop correction 
is small, while the correction in each decay width is large. With the experimental value 
Tfe = 37.04 we can estimate as{mb) as 

as{mb) = 0.18. (34) 

which is close to the value of as by running down from the scale /i = Mz- However, the 
relativistic correction is not significant in the determination of a^. Similar estimation can 
be done for J/\l/, the determined asifric = 1.6GeV) is 0.195. It is smaller than the expected. 
This indicates that effects in higher orders of v may be still significant. It also should be 
noted that the pole mass rric is not known precisely as mi). 

To summarize: In this work we tried to define the difference between a quarkonium 
mass and the twice of the heavy quark mass in terms of NRQCD matrix elements. We 
separated the effect introduced by the difference in relativistic corrections for several decays 
of quarkonium and find that the relativistic correction is determined by the difference. With 
the determination decay widths considered here are enhanced by the relativistic correction. 
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